Nonorthogonal pairs of copropagating optical modes in deformed microdisk cavities 
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Recently, it has been shown that spiral-shaped microdisk cavities support highly nonorthogonal pairs of co- 
propagating modes with a preferred sense of rotation (spatial chirality) [Wiersig et ai, Phys. Rev. A 78, 053809 
(2008)]. Here, we provide numerical evidence which indicates that such pairs are a common feature of de- 
formed microdisk cavities which lack mirror symmetries. In particular, we demonstrate that discontinuities of 
the cavity boundary such as the notch in the spiral cavity are not needed. We find a quantitative relation between 
the nonorthogonality and the chirality of the modes which agrees well with the predictions from an effective 
non-Hermitian Hamiltonian. A comparison to ray-tracing simulations is given. 

PACS numbers: 42.25.-p, 42.55.Sa, 05.45.Mt, 42.60.Da 



I. INTRODUCTION 

Optical microcavities allow trapping of photons for a long 
time Tc in very small volumes This enables one to 

control light and matter at the nano- and microscale, which 
is important for many applications such as ultralow thresh- 
old lasers iS [H and single-photon emitters |0, IH. Mi- 
crodisks 10,1711, microspheres and microtoroids iflollTTll 
support so-called whispering-gallery modes. Due to total in- 
ternal reflection of the photons at the boundary of the cavity 
these modes have very high quality factors Q = ujTc, where 
u! is the resonance frequency. However, as a consequence of 
the rotational symmetry of a microdisk the in-plane light emis- 
sion from these modes is isotropic, which is a disadvantage for 
many applications. This problem can be solved by deforming 
the boundary of the cavity lfl2l - [l4|] . Several shapes for unidi- 
rectional light emission have been proposed, for example the 
spiral cavity Iil5i - tl7i] . the annular cavity [18], the limagon cav- 
ity lfl9ll . the circular disk with a point scatterer ||20| . and the 
notched ellipse fHI . 

The spiral cavity is a well-studied system ifTsl - flTll . In 
polar coordinates the boundary of this cavity is defined as 
= R{1 — ^(f) with deformation parameter e > and 
radius R at (f> = 0. The radius jumps back to i? at = 2tt, 
creating a notch. The spiral cavity appears to be special in 
the list of studied geometries for two reasons: (i) it lacks any 
discrete spatial symmetry and (ii) the boundary curve exhibits 
a singularity. In Refs. ll22il23ll it has been demonstrated that 
the modes in this open system come in highly nonorthogonal 
pairs. Moreover, each pair of modes shows a strong spatial 
chirality, in the sense that both modes have mainly counter- 
clockwise (CCW) propagating components, while the clock- 
wise (CW) component is weak in both modes. (It is im- 
portant to emphasize that our usage of the term "chirality" 



should not be confused with optical activity in chiral media, 
see, e.g., ll24ll .) The appearance of nonorthogonal and chiral 
modes in the spiral cavity has been traced back to the asym- 
metric scattering between CW and CCW propagating waves 
at the notch |22]. 

The aim of the present paper is to show that this effect 
also appears in cavities without boundary singularities. We 
consider two different geometries. One is a rather represen- 
tative example, an asymmetric version of the limafon. The 
other example is a variant of a curve of constant width ||25| . 
This exotic geometry helps to clarify the role of ray dynamics. 
Our results indicate that the nonorthogonal and chiral pairs of 
modes appear in any cavity geometry which is a sufficiently 
small deformation of the circle and lacks any mirror symme- 
tries. It should be mentioned that in experiments fabrication 
tolerances introduce small asymmetries quite naturally, see, 
e-g-, lUS- In this sense, the asymmetric shapes are generic 
which is, however, not reflected in the list of geometries stud- 
ied in the literature. 

This paper is organized as follows. Section reports our 
numerical results on the properties of optical modes in the 
asymmetric limagon cavity. In Sec.|lll]we introduce an effec- 
tive non-Hermitian Hamiltonian which describes the relation 
between nonorthogonality and chirality. The ray dynamics 
in the asymmetric limagon cavity is presented in Sec. |IV] In 
Sec. [V] we discuss the cavity with a boundary curve of con- 
stant width. We summarize our results in Sec. |VT] 



II. MODES IN THE ASYMMETRIC LIMA^ON 

In the case of (deformed) microdisk cavities with a piece- 
wise constant effective index of refraction n{x, y). Maxwell's 
equations can be reduced to a two-dimensional scalar mode 
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where w = cfc is the frequency, k is the wave number in 
vacuum (outside the microdisk), and c is the speed of hght 
in vacuum. The mode equation ([T]i is vaHd for both trans- 
verse magnetic (TM) and transverse electric (TE) polariza- 
tion. For TM polarization the electric field E{x,y,t) cx 
(0, 0,Re[ip{x, i/)e~*'^*]) is perpendicular to the cavity plane. 
The wave function t/j and its normal derivative are con- 
tinuous across the boundary of the cavity. For TE polariza- 
tion, ip represents the z-component of the magnetic field vec- 
tor Hz- Here, the wave function il' and n{x, y)^'^duijj are con- 
tinuous across the boundaries [|27ll. At infinity, outgoing wave 
conditions are imposed, which results in quasibound states 
with complex frequencies lj in the lower half-plane. The real 
part is the usual frequency and the imaginary part is related 
to the lifetime Tc = — l/[2Imc<;] and to the quality factor 
Q = -Recj/plmw]. 

In polar coordinates the boundary shape studied in this sec- 
tion is given by 



p{<p) = R[l + El cos 4> + 62 C0S(2(/) + S)] 



(2) 



with (5 e [0, 2tt). This is illustrated in Fi g.fT] The special case 
£2 = is the limagon cavity lfl9ll26il28l - j33ll . For nonzero ei, 
£2 and (5 7^ 0, TT the system does not possess any mirror sym- 
metry. We therefore call the cavity described by Eq. (|2) the 
asymmetric Umagon. The ray dynamics in such a geometry is 
mainly chaotic as in the case of the symmetric limagon. The 
asymmetric version can be considered as a typical smooth de- 
formation of a circular disk, as the sum of the first three terms 
of a Fourier expansion of an arbitrary periodic function can be 
written as in Eq. (|2|l. 

If not stated otherwise, in the following we use 5 = 
sa 0.6187r, ei = 0.1, and 62 ^ 0.075. We have cho- 
sen the golden ratio for Sj-Kio ensure that we are not too close 
to the symmetric situations (5 = and tt. The effective index 
of refraction is set to n = 3.3 (e.g. GaAs) and we only con- 
sider TM polarization. To compute the complex frequencies 
and wave functions of quasibound states we use the boundary 
element method Il34ll a nd, for comparison, the wave matching 
method, see, e.g., OSll . We always find very good agreement 
between these two methods. 

Figure |2] shows the resonances in a typical region of the 
complex plane of normalized frequencies Q, = loR/c = kR. 
It can be clearly seen that modes here always appear in nearly 
degenerate pairs, even though the deformation of the cavity 
boundary is not that small. This is due to the rather weak 
coupling of CW and CCW propagating waves in these open 
disk-like cavities. One could guess that one member of such 
a pair is a CW-propagating mode and the other one is a CCW- 
propagating mode. But this is not true, as we will see in the 
following. 

Figure [3] depicts a typical example of such a pair of nearly 
degenerate modes. The splitting in real part ARe(51) w 
4.4 X 10~^ and imaginary part Alm(fi) w 10~'^ is very small. 
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FIG. 1: (Color online) The left panel shows the boundary 
parametrization p((l>) in Eq. ^ with ei = 0.1 and £2 = 0.075, 
for different values of 5. The curves are shifted vertically for bet- 
ter comparison. The right panel shows a top view of the cavity with 
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FIG. 2: Position of dimensionless complex resonance frequencies 
with TM polarization for the cavity in Eq. (|2) with refractive index 
71 = 3.3 and shape parameters £1 = 0.1, £2 = 0.075, and 5 = 
Open circles (crosses) mark the slightly higher-Q (lower- 



Q) mode of a given pair of modes. 



The spatial mode pattern is difficult to distinguish by eye. A 
closer look at the far-field pattern shown in Fig. 2]reveals that 
they have the same envelope, but there clearly are different 
oscillations on top of this envelope. Note that not only the far- 
field pattern indicates unidirectional light emission, the near- 
field pattern in Figs.|3lc) and|3jd) shows that the emission fol- 
lows a single-lobe beam, which is interesting for applications. 



Following Refs. 11221 12311 we analyze the mode pattern by 
expanding the wave function inside the cavity in cylindrical 
harmonics. 
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where J,„ is the mth order Bessel function of the first kind. 
Positive (negative) values of the angular momentum index m 
correspond to CCW (CW) traveling-wave components. As the 
origin of this expansion we choose the center of mass of the 
cavity, {x,y) = {eiR/2,0). Note that this particular choice 
does not affect our conclusions (as long as the origin is cho- 
sen inside the cavity). The coefficients in the expansion in 
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= {am'' ± am' ) /2 being CW and CCW traveling-waves, 
respectively, as can be seen in Fig. |5jd). However, the CCW 
superposition has a much larger amplitude. It is important to 
emphasize that these superpositions are not eigenmodes of the 
cavity as they are composed of two modes with slightly dif- 
ferent frequencies and Q-factors. 
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FIG. 3: (Color online) Intensity \tpf of the nearly degenerate pair 
of modes in the asymmetric limagon with (a) Qi = 12.319807 — 
iO. 00089 and (b) ^2 = 12.319851 - i0.0009; cf. Fig.[2] The same 
color map has been used in panels (a) and (b). Panels (c) and (d) 
shows the corresponding exterior mode pattern at some small dis- 
tance away from the cavity (white region). The far-field patterns are 
shown in Fig.|4] 
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FIG. 4: Far-field patterns of the modes in Fig. [5] The solid and 
dashed lines coiTespond to the modes in Figs.|3}a) and|3lb), respec- 
tively. For the definition of the far-field angle see Fig.[T] 



Eq. (|3]l are naturally given in the case of the wave matc hing 
method Issll . In the case of the boundary element method ll34ll 
we use a Fourier transformation of the wave function to deter- 
mine the coefficients. Both approaches give identical numeri- 
cal results. 

In Fig. 13 a) we can observe that for both modes the angu- 
lar momentum distribution |q;„iP is dominated by the CCW 
component, i.e., none of the two modes can be classified as a 
CW traveling-wave mode. The small difference between the 
expansion coefficients of the modes can be seen in Figs. 13b) 
and|5]!c). For negative angular momentum index both the real 
and the imaginary part of am have a different sign for the two 
modes. That means that we can construct superpositions with 
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FIG. 5: (Color online) Angular momentum distributions (black 
solid line) and ain (green dashed) of the modes in Fig. |3] (real part 
normalized to 1 at maximum): (a) absolute value squared, (b) real 
and (c) imaginary part, (d) superpositions 



(black solid) and a„ = {a 
a factor of 6). 
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It is convenient to use the angular momentum representa- 
tion (O to define the (spatial) chirality of a mode by 



max 



\ Z^m— 



2 y= 

1 Zl^r 



(4) 



If the weight of the CW and CCW components is equally 
distributed then the chirality is a = 0. This is for instance 
the case for a cavity which possesses a mirror symmetry. To 
see this, choose the coordinate system such that p(— 0) = 
In the angular momentum representation (|3) modes with 
positive (negative) parity —0) = ±^(/3, 0) must have 
a_,„(— 1)™ = iflm- In both cases, lo-mp = IfljnP which 
according to Eq. (|4]i gives a chirality a = 0. The same is true 
in the case of a closed system (in nonlinear dynamics a closed 
cavity is called a billiard) with real-valued frequencies, where 
a'*ij„(— 1)'" = am which again leads to |a_mP = |amP- 
Note that the former statement is correct for the special case 
of the circular microcavity with its degenerate pairs of modes 
only if the linear superpositions leading to standing waves are 
chosen. The other extreme of full chirality, a = 1, is realized 
if a mode has no CW (or CCW) component at all. 

The asymmetric limafon (with 8 ^ 0, tt) lacks any mirror 
symmetry. Moreover, the system is open and has complex- 
valued frequencies due to the outgoing-wave conditions at in- 
finity. For the modes in Figs. |3la) and[3jb) we find numeri- 
cally a Ri 0.839 and a w 0.8404, respectively. Hence, both 
modes show a strong chirality. 

Now we demonstrate that the modes not only have a strong 
chirality but that they are also pairwise highly nonorthogonal. 
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FIG. 6: Time evolution of CCW (solid lines) and CW (dashed) 
components in semilogarithmic scale. The traveling waves are su- 
perpositions of the modes with frequencies Qi and SI2 depicted in 
Fig. [3] The upper (lower) panel contains the dynamics starting with 
a pure CCW (CW) traveling wave. Time is measured in units of 
Ti = 27r/Reni. 



To quantify the nonorthogonality we compute the normalized 
overlap integral of two modes ipi and ^2 over the interior of 
the cavity C (see also, e.g., ll36ll ) 



S 



I L dxdy ijjl-ip2 



^ dxdy ■ipl'ipi^ dxdy ■4>*'ip2 



(5) 



In the case of orthogonal states S* = and in the case of 
collinear states S = 1. It is easy to show that in the pres- 
ence of a mirror symmetry or for a closed system the over- 
lap S vanishes. For the pair of modes in Fig. |3] however, we 
find S w 0.7236 reflecting a strong nonorthogonality. Mode 
nonorthogonality is important as it implies excess quantum 
noise 1 37-41]. 

One might think that the strong chirality must have a large 
influence on the dynamics of waves in such a cavity. For the 
asymmetric lima^on the impact is, however, weak. The upper 
panel of Fig.|6]shows the dynamics of the CCW and CW com- 
ponents of an initially (at time t = 0) CCW traveling wave 
using proper superpositions of the modes shown in Fig. |3](we 
follow the procedure explained in detail in Ref. ||22| ). Only a 
very weak scattering into the CW component can be observed. 
The lower panel of Fig. |6]shows the situation where the initial 
wave is purely CW propagating. Here, the backscatteiing into 
the CCW component is significantly larger, though still weak. 

The chirality of the modes in Fig. [3] as a function of the 
asymmetry parameter 6 in the interval [0, it] is shown in the 
lower panel of Fig. [T] The curves have been computed by 



starting with S = n 



V5-1 



decreasing S in small steps, and 



thereby following the modes. Second, S has been varied be- 
tween TT ^2^^ and TT. 200 discretization points on the (5-axis 
have been used for each mode. As the total change in the fre- 
quency is rather small in this range of parameter variation we 
do not observe any avoided resonance crossing. For 5 = and 
6 = TT the system possess a mirror symmetry and therefore we 



observe no chirality, a = 0. The maximum chirality of about 
0.845 is attained at (5 « 2, which is close to the value that we 
mainly use in this paper, (5 = tt w 1.94. The upperpanel 
shows the level splitting Ail = |Re(rii) — Re(ri2)| and the in- 
dividual decay rates of the two modes — Im(f2i) as a function 
of the asymmetry parameter 5. Note that the level splitting is 
always much smaller than the individual decay rates, i.e., the 
spectral width of the two modes strongly overlap. 
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FIG. 7: Upper panel: individual decay rates — Im(Sli) (solid and 
dashed line) and the level splitting Af2 = |Re(f2i)— Re(5l2)| (dotted 
line, scaled by a factor of 10) of the pair of modes in Fig.|3]vs shape 
parameter S. Lower panel: corresponding chirality a as function 
of 5. Note that two curves are on top of each other. 

Another example of a pair of nearly degenerate modes is 
shown in Figs.[8]and|9l Again we find that both modes exhibit 
a strong chirality (a w 0.8793 and a « 0.8658) and a signif- 
icant nonorthogonality [S « 0.7778). This is in particular re- 
markable as the quality factor of the modes is about 3.5 x 10^. 
For such enormously high quality factors one would expect a 
behavior similar to that of orthogonal states in a closed sys- 
tem such as in a billiard. This reasoning is, however, too 
naive since the frequency splitting is much less than the in- 
dividual decay rates. The resonances therefore strongly over- 
lap, which is usually considered as a feature associated with 
strongly open systems. 

It is also important to mention that the overlap of long-lived 
modes from different pairs is significantly smaller For in- 
stance, for one mode in Fig.|3]and one from Fig.|8]we always 
find 5* < 6 X 10^^. That means the nonorthogonality is sig- 
nificant only within each pair of long-lived modes. 

Figure \W\ summarizes the results on the chirality and the 
pairwise nonorthogonality of the modes which are present in 
the frequency regime considered in Fig.|2] The mode pair with 
the highest Q-factors (fli « 12.0960925 - ilM x 10""^ and 
n2 ~ 12.09609251 - il.ll x 10"'^) is, however, not shown 
here. In this (and only in this) case the numerical computation 
of the chirality is not fully converged due to the exceptionally 
strong degeneracy of this particular mode pair even though we 
use up to 32 000 discretization points in the boundary element 
method. From Fig. [TO]it can be observed that the chirality 
and the overlap are correlated. This correlation can be ex- 
plained by an effective non-Hermitian Hamiltonian, which is 
discussed in the next section. Finally, we remark that 29 of 
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FIG. 8: (Color online) Intensity of the nearly degenerate pair 
of modes in the asymmetric lima§on with (a) f2i = 12.73070292 — 
il.83 X 10"^ and (b) 0.2 = 12.73070286 - il.88 x 10"^ ; cf. Fig.|2] 
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FIG. 9: (Color online) Angular momentum distributions Ofm' (black 
solid line) and Ofm' (green dashed) of the modes in Fig.[8]normalized 
to 1 at maximum: (a) absolute value squared, (b) real and (c) imag- 



inaiy parts, and (d) superpositions q+ — (om' + (black 



solid) and a„ 
of 8). 



(a 



(1) 



J2) 



)/2 (red dashed, multiplied by a factor 



the 3 1 considered pairs of modes exhibit a larger CCW com- 
ponent. Only two pairs have a larger CW component. 

For smaller frequencies the chirality and the nonorthogo- 
nality are weaker (not shown). An intuitive explanation is that 
for very small frequencies, i.e., large wavelengths, the modes 
do not feel the asymmetry of the boundary shape anymore. 



III. EFFECTIVE NON-HERMITIAN HAMILTONIAN 



Reference 112211 introduced a simple toy model to describe 
the main features of the chirality and nonorthogonality of 
modes in the spiral cavity. Here we use the same two-by-two 
non-Hermitian and nonsymmetric matrix 



H 



LOq 

OJQ 



r 

'qV* 



V 

r 



(6) 



For the convenience of the discussion, our interpretation is 
that this matrix describes the dynamics of the wave function ij.' 
in slowly-varying envelope approximation ll42ll in the time do- 
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FIG. 10: Chirality a vs spatial overlap S of pairs of almost degen- 
erate modes in the asymmetric limafon. Open circles (crosses) mark 
the slightly higher-Q (lower-Q) mode of a given pair of modes com- 
puted numerically from Maxwell's equations; cf. Fig. |2] The solid 
line is the analytical prediction of the theoretical model, Eq. il5\ . 



main by a Schrodinger-type equation 
. d 



(7) 



When deriving this equation from the Maxwell's equations 
one assumes that the optical field varies slowly in time (not 
necessarily in space) with respect to a reference frequency 
which we choose to be close to the two nearly degenerate 
modes of interest. 

The nonsymmetric Hamiltonian matrix (|6]l is defined in the 
CCW/CW traveling-wave basis 



t2 = 



(8) 



The eigenvectors of the first matrix in Eq. (|6) on the right-hand 
side belong to the CCW and CW traveling waves with, for 
simplicity, equal frequency G C in the absence of any cou- 
pling between them. The second matrix accounts for coupling 
of CCW and CW traveling components. The diagonal ele- 
ments are given by the total decay rates and frequency shifts 
r e C which are assumed to be equal for simplicity. The off- 
diagonal element V = |V^|e''^ G C describes scattering from 
a CW traveling wave into the CCW traveling wave. The other 
off-diagonal element 77^* describes scattering from a CCW 
traveling wave into the CW traveling wave. The latter scatter- 
ing is assumed to be weaker, i.e., \ri\ < 1. Therefore, here \ri\ 
plays the role of the asymmetry parameter 

Note that a standing wave basis can be chosen as 



1 
71 



S2 



V2 



-1 
1 



(9) 



While the traveling-wave basis dHJ corresponds to terms 
with rn > and 7ti < in the angular momentum representa- 
tion ([3]), the standing wave basis corresponds to cos (mcj)) and 
sin (mcj)) with nonnegative m. In the standing wave basis the 



6 



non-Hermitian Hamiltonian matrix is symmetric (as required 
by time-reversal symmetry): 



H = 



u!o + r - 



V + 7jV 



V+rjV' 



(10) 



The complex eigenvalues of the matrix (|6]l are given by 

io± = Loo + r±^j\v\- (11) 

The (not normalized) right-hand eigenvectors in the 
CCW/CW traveling-wave basis turn out to be 



a± 



1 



±076 



(12) 



These eigenvectors explain the mode structure observed in 
Sec. Hn including the sign difference in the CW components 
as well as the relative weight of the CCW and CW compo- 
nents. The weight of the first component (corresponding to 
CCW traveling waves) squared ~ 1 is much larger than that of 
the second component (corresponding to CW traveling waves) 
squared ^ \r]\, cf. Figs.|5]and|9l Hence, the 2x2 model pre- 
dicts an identical chirality for both modes, 



1 



1^1 



(13) 



which is nonzero in the case of asymmetric coupling ( I ?7 1 7^ 1). 
Note that the chirality a does not depend on the coupling 
strength \V\. As a consequence, even an infinitesimal cou- 
pling can lead to a significant chirality. The reason behind 
this singular behavior is the two-fold degeneracy of the un- 
perturbed modes. 

The two eigenvectors in Eq. (fT2] i are, in general, nonorthog- 
onal, i.e., the normalized overlap 
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1 



\v\ 



1 



\v\ 



does not vanish. Using this result and Eq. ( fT3] l we arrive at 
a relation between the chirality of the two modes and their 
overlap. 



2S 

1 + S 



(15) 



Figure [TO] compares this prediction with the data obtained 
from numerical solutions of Maxwell's equations. It can be 
seen that the 2x2 model works very well. 

For the case ?; = 0, the Hamiltonian (|6]l exhibits an ex- 
ceptional point ll38[|43ll . i.e., not only the eigenvalues become 
degenerate. 



iu± = ojQ +T 



(16) 



but also the eigenvectors collapse to a single one, which in the 
CCW/CW traveling-wave basis reads 



a± = 



(17) 




lm(r|) 



Re(ii) 



FIG. 1 1 : (Color online) Complex-square-root topology with a branch 
point singularity at the exceptional point of the Hamiltonian 



The corresponding complex-square-root topology of the 
eigenvalues ( fTTI ) with a branch point singularity is shown in 
Fig.[TT] The eigenvector at the exceptional point in the stand- 
ing wave basis ^ is given by 



a± 



+ i 



(18) 



This is a chiral state in the sense ofRefs. Hlil. 

The solutions (V'i(t), ip2{t)) of the Schrodinger-type equa- 
tion (|7| can be found analytically. We consider the time scale 
related to the strength of the mode splitting 



T 



1 



^^-j and the decay time 



1 



lm(a;o + T) 



(19) 



(20) 



In the following we restrict ourselves to the case r <C T and 
t <C T which is the relevant regime for small \t]\. We find for 
a wave propagating initially at t = in CCW direction with 
normalized amplitude 



iV'il 



't/r 



(21) 
(22) 



The component IV'il decays from 1 to in an exponential 
manner and \ip2 \ increases from zero to the value 



1^2 



T 



(23) 



and then decays to zero. For I77I ^ 1 the value of | '021 max 
is very small. This corresponds to the situation in the upper 
panel of Fig.|6l Note that Eqs. (I2ni-(|22]| are exact for 77 = 6, 
i.e., at the exceptional point; see also 1221 14611 . 
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For a wave propagating initially at t = in CW direction it 
follows that 



1^1 



1^2 1 = e-*/^ 



1 t 



T 



-t/T 



(24) 
(25) 



The component \tp2\ decays from 1 to in an exponential 
manner However, \ increases from zero to the value 



(26) 



and then decays to zero. Here, |V'i|max, in principle, can be 
close to 1, as |r/| is a small number The ratio 



l^2|max |?7| 



(27) 



reflects the fact that the scattering from CW to CCW is l/\ri\ 
times stronger than the scattering from CCW to CW. 

For the modes in the asymmetric limagon in Fig. |3] we es- 
timate t/T w 0.024, from Fig. |5] I77I w 1/6, and there- 
fore IV'ilmax ~ 0.022 and |V'2|max ~ 0.0037, in reasonable 
agreement with Fig. |6] For the modes in Fig. |8]we estimate 
iV'ilmax ~ 0.017 and |V'2|max « 0.0021. If we apply the 
same analysis to the spiral cavity then we find that the | V'l |max- 
values are significantly larger This is consistent with the fact 
that the mode splitting and the chirality in the spiral cavity is 
larger. 



IV. RAY DYNAMICS IN THE ASYMMETRIC LIMAgON 

Figure [12] shows the ray analog of Fig. |6] In the up- 
per (lower) panel a bunch of 12 500 rays in the asymmetric 
limagon propagating in CCW (CW) direction with angle of 
incidence |x| distributed uniformly well above the critical an- 
gle for total internal reflection, | sin(xc)| = 1/n w 0.3, has 
been launched. In our case, we select | sinx| > 0.5 in or- 
der to eliminate the very short-lived rays. It can be seen that 
the scattering from CCW to CW is weaker than from CW to 
CCW, as in the wave calculations in Fig.|6l However, the ratio 
of the scattering rates, corresponding to l/lr/j, is below 2, i.e., 
the asymmetry of scattering seems to be weaker for the ray 
dynamics. Note, however, that the asymmetry here depends 
on the initial conditions of the rays. Restricting the initial rays 
to regions near confined periodic ray trajectories can enhance 
l/\ri\ to around 10. 

Having observed this correspondence of rays and waves in 
terms of scattering, we show now that there is no such corre- 
spondence of rays and optical modes in terms of spatial chi- 
rality. To see this, consult Fig. [13] which is the ray analog 
of Fig. 13 a). In total 40 000 rays with initially uniform dis- 
tribution along the boundary of the cavity and uniformly dis- 
tributed |sinx| G (0.5,1) have been started. The quantity 
sinx is related to the angular momentum of a ray. In the 
case of the circle the relation is given by nkRsinx = m 
with angular momentum index m ITsll . Positive values of 




■3 1 



3 0.1 r 



0.01 







CW 






CCW ^ 






1,1,- 



100 200 300 400 

length of ray trajectory in units of R 



500 



FIG. 12: Time evolution of amplitude (defined as square root of in- 
tensity) in semilogarithmic scale coixesponding to CCW (solid lines) 
and CW (dashed) propagating light rays in the asymmetric limafon. 
The upper (lower) panel shows the dynamics starting with a set of 
pure CCW (CW) propagating rays in analogy to the wave dynamical 
considerations in Fig.|6] Time is proportional to the geometric length 
of ray trajectories. 



sin X correspond to CCW propagation direction and negative 
to CW propagation direction. After time t = 50, measured 
in length of ray trajectory in units of R, the remaining inten- 
sity of rays approaches a survival probability distribution ||47|l. 
This distribution reflects the long-time behavior of the light 
rays which can be compared to the properties of the long- 
lived modes. The survival probability distribution as function 
of sin X is plotted in Fig. [T3] We do not observe a chirality in 
this distribution, i.e., the amount of CW propagating rays is 
with 49.8 per cent roughly equal the amount of CCW propa- 
gating rays, which is in strong contrast to the properties of the 
optical modes in Fig. Oa). 



CW 



CCW 



-1 -0.5 0.5 1 

sin(x) 

FIG. 13: Intensity of long-lived rays vs angle of incidence x in 
the asymmetric limajon. Positive (negative) x correspond to CCW 
(CW) propagation direction. 

Neither does the spatial chirality show up in an extended 
ray dynamics including first-order wave corrections such as 
the Goos-Hanchen shift (GHS) and the Fresnel filtering (FF). 
The GHS is a lateral shift of totally reflected beams along the 
optical interface [48], i.e., the points of incidence and reflec- 
tion do not coincide. In the case of the FF lH^-jUt], partial 
waves with angles of incidence below the critical angle for 
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total internal reflection are (partially) refracted out of the cav- 
ity, leading to a shift Ax of the partial waves between the 
incident and outgoing angles - i.e., a violation of Snell's law. 
In the short-wavelength limit A the GHS and FF disap- 
pear leading to the standard ray dynamics of geometric optics. 
Such wave corrections have been used to explain properties 
of optical modes in deformed subwavelength-scale microdisk 
cavities [26, 52, 53j|. 

The presented ray simulations for the asymmetric limagon 
show that the weak asymmetry in the scattering between CCW 
to CW propagating rays does not lead to a chirality in the sur- 
vival probability distribution. Hence, the ray dynamics cannot 
explain the chirality observed in the optical modes. To sup- 
port this finding, in the next section we discuss another cavity 
geometry, for which it is rigorously proven that there is no 
scattering between CCW and CW propagating rays and there- 
fore no asymmetry in the scattering of rays. Nevertheless, the 
modes show significant chirality and nonorthogonality. 

V. GUTKIN'S BILLIARD OF CONSTANT WIDTH 

Gutkin studied a class of convex billiards (closed cavities) 
of constant width ll25ll . i.e., for any point at the boundary the 
maximal distance to other points of the boundary is a con- 
stant; see, e.g.. Fig. [141 The (conventional and extended) 
ray dynamics in such billiards is characterized by a phase 
space which is strictly separated into two parts correspond- 
ing to CW and CCW motion, i.e., there is no scattering from 
CW to CCW propagating rays. Optical microcavities of con- 
stant width have been studied in the context of directional light 
emission ifsill . The parametrization of the class of boundary 
shapes in the (x, ?/)-plane is most conveniently given in the 
complex variable z ^ x + iy 

z{a) = z(0) - ^ ^ (e»("+i) - l) (28) 

with a £ [0, 2tt), a-n = ai = 0, and a2n = for n > 0. 
Here we consider a realization without mirror symmetry and 
where the CW and CCW component in phase space are almost 
fully chaotic: z(0) = (1/4 - i)R, ao = R, = iR/8, = 
(1 + i)R/4, and a2k+i = for k > 2. The constant width is 
W ~ 2R. This boundary curve is illustrated in Fig. [14] 

We find that the optical modes in a microcavity with such a 
boundary shape also appear in nonorthogonal and chiral pairs 
of modes. Figures[T5]and[T6]show the mode pattern and its an- 
gular momentum decomposition (origin is the center of mass) 
of a typical pair of modes for TM polarization and refractive 
index n = 3.3. The chirality turns out to be a w 0.726 and 
w 0.717. The spatial overlap of both modes is around 0.56. 
These values are in good agreement with the result from the 
effective Hamiltonian in Eq. JTSl l. 

A proper rotation of Fig. [T4]shows that the boundary curve 
is less asymmetric than the asymmetric limagon. We there- 
fore consistently observe smaller chirality a and overlap S for 
the cavity geometry defined by Eq. ( i28] l. Moreover, we find 
that the relative number of CW copropagating pairs is larger 
if compared to the asymmetric limagon. 








X 

FIG. 14: Billiard boundary of constant width as defined in Eq. J28t 
with z{0) = (1/4 - i)R, ao = i?, as = iR/8, as = (1 + i)R/4, 
and a2k+i = for fc > 2. The straight lines show four possible 
ways (out of infinitely many) to measure the width W. The value of 
W is always 2R. 




FIG. 15: (Color online) Intensity \^p\'^ of the nearly degenerate pair 
of modes in Gutkin's billiard of constant width with (a) Qi — 
12.386847 - iO. 00561 and (b) ft2 = 12.386995 - iO. 005695. 



In a billiard of constant width there is no scattering from 
CW and CCW propagating rays and vice versa. Using this 
fact, we can rule out the asymmetric scattering of rays as the 
origin of the spatial chirality in this system. 



VI. SUMMARY 

The nonorthogonality and spatial chirality of mode pairs in 
two asymmetrically deformed microdisk cavities, the asym- 
metric limagon and Gutkin's cavity of constant width, has 
been studied. Our results indicate that the appearance of such 
nonorthogonal chiral pairs is a common feature of deformed 
microdisks which lack mirror symmetries. Using an effec- 
tive non-Hermitian Hamiltonian we have linked these two in- 
teresting effects and explained them by the asymmetric scat- 
tering between clockwise and counterclockwise propagating 
waves, expressed by an asymmetry parameter \ri\. We have 
shown that the observation of these effects in dynamical ex- 
periments with waves depends not only on the ratio between 
the period related to the mode splitting T and the decay time t 
but also strongly on the asymmetry parameter |77|. Finally, we 
have demonstrated for the considered cavities that there is no 
significant chirality in the survival probability distribution of 
rays. This is in strong contrast to the case of the spiral cavity. 
This observation shows that the nonorthogonality and spatial 
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FIG. 16: (Color online) Angular momentum distributions Qm' (black 
solid line) and (green dashed) of the modes in Fig.|15|normal- 
ized to 1 at maximum: (a) absolute value squared, (b) real and (c) 
imaginary parts, and (d) superpositions — (ofm' +Qm')/2 (black 



solid) and a„ 
of 3.5). 



(a 



(1) 



„(2) 



)/2 (red dashed, multiplied by a factor 



chirality is in general a wave dynamical effect. 



It remains an interesting question for future research to ask 
for the size of the chiral effects in the case of circular-shaped 
disks with surface roughness arising from the inevitable im- 
perfections in the fabrication process. In such a case we also 
expect the appearance of nonorthogonal pairs of copropagat- 
ing modes. But the nonorthogonality and the chirality might 
be small for the case of weak surface roughness that can be 
achieved nowadays in state-of-the-art experiments. Moreover, 
averaged over many mode pairs the chirality should be close 
to zero due to the random character of the boundary profile. 

We believe that our results are not only important for de- 
formed microdisks but also for other types of optical micro- 
cavities (microspheres and microtoroids) and for open quan- 
tum (wave) systems in general. 
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